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اندر کنش سد و مخزن در قلمرو زمان 
تاد هی ۱ رضا عطارنةا(؟ 


چکیده در اين مقاله تحلیل دقیق پویای اند رکنش سد و محزن در قلمرو زمان انجام می‌شود. ضخامت بدنه سد. متغیر و سازه سد, انعطاف پذیر و 
محزن تا بی‌نهایت ادامه دارد. شرایط مرزی تشعشعی در مخزن و تغییر شکل بدنه سد به دقت در رابطه‌ساز ی گنجانده شده‌اند تا همه بعدهای 
پدیده فیزیکی اند رکنش سد و مخزن را دربربگیرند. فشار هیدروپویای در پشت سد ازمعادله فشار موج پیروی می‌کند. الگوسازی سازه سد به 
صورت اویلر- برنولی و با فرض قائم بودن بدنه در بالادست انجام پذیرد. روش نوینی برای تحلیل تیر با مفطع متغیر ارایه می‌شود. این روش وین 
بر مبنای کاربرد تابع‌های نوینی است که تابم‌های شکل مبنا( 01۳ ) نام دارند. این تابع‌ها از حل معادله دیفرانسیل حرکت تیر اویل ر- برنولی با 
مقطع متغیر به دست م یآیند و با کمک آن‌ها تابع‌های شکل پویای نوینی حاصل می‌شوند که در تحلیل تیرهای با مقطع متغیر به کار می‌روندد. 
پاسخ پویای سد و محزن در برابر زلزله ال سنترو بررسی و تغییر مکان سد در حالت محزن پر وحالی محاسبه شده است. مقایسه ننیجه‌های حاصل 


با نتیجه‌های موجود در نوشته‌های فنی» نشان‌دهنده شایستگی و دقت روش ارایه شده است. 
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مقد مه 

برای بررسی پاسخ پویای سد در برابر زلزله بایستی 
سامانه سد و مخزن را هم‌زمان مورد تحلیل و بررسی قرار 
داد انر ماه شتا سل ما زمرت و هاش تسس 
می‌باشد. روش‌های گوناگونی برای تحلیل سازه سد و 
مخزن به کار می‌روند. در ایین بخش. نخست درباره 
بو هويم کار گنس یی انا 
توضیح داده می‌شود و سپس روش نوینی برای تحلیل 
اندرکنش سد و مخزن معرفی خواهد شد. 

توزیع فشار هیدرودینامیک را برای نخستین بار 
وسترگارد در سال ۱۹۳۳ بر روی سدهای صلب ارایه داد 
[1] در سال ۸۹7۷ 0۳2۵ یک روش تحلیلی برای 
توزیع فشار هیدرودینامیک روی یک سد صلب با بدنه 
قائم پيشنهاد کرد [2]. در سال ۱۹۷۸ 0۳۵08 با 
چشم‌پوشی از فشردگی پذیری سیال, رابطه‌سازی فشار 
هیدرودینامیک برای سدهای صلب با دیواره بالادست 
مایل با زاویه ثابت را ارایه داد [3]. 1۸6 در سال ۱۹۷۹ 
یک روش دقیق برای حل مسأله اندرکنش سازه و سیال 
در قلمرو بسامد به‌دست آورد [4]. در سال ۱۹۸۲ ]با 
گسترش رابطه‌سازی ۳۷۵08 یک روش دقیق برای 
سدهای صلب و بدنه شیب‌دار و مخزن مثلثی ارایه 
کرد[5]: 

تحلیل اندرکنش سد و مخزن با روش‌های عددی 
مختلفی انجام می‌پذیرد. روش‌هایی مانند جزء‌ی محدود. 
جزء‌ی مرزی» جزء‌ی نامحدود و روش تفاضل محدود از 
جمله روش‌هایی هستند که در تحلیل مسأله اندرکنش 
مورد توجه قرار گرفته‌اند. ۳۵۳۷65 و 00072 در سال 
۶ یک برنامه رایانه‌ای بر مبنای روش جزءی محدود 
برای تحلیل مسأله اندرکنش سد و مخزن ارایه دادند [6]. 
در سال ۱۹۹۰ عحعط) و ۷1۲ تحلیل لرزه‌ای سد- مخزن 
و فونداسیون را با روش جزءی مرزی انجام دادند[7]. در 
سال ۱۹۹۱ 50047 و ۸۵۸65 با ترکیسب روش جزءی 
محدود و جزء‌ی مرزی یک روش بر اساس گام‌های 
زمانی برای تحلیل اندرکنش سد - مخزن پيشنهاد کردند 
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[8]. ۲۵۲۱ و 1600 در سال ۲۰۰۳ برای تحلیل دو 
بعدی آندرکنش سد- مخزن» روشی بر مبنای جزء‌ی 
محدود و نامحدود پایه نهادند [9]. در سال 
۶ 16100012792120 انْ_درکنش سسد- مضزن- 
فونداسیون را در قلمرو زمان با روش تقابل دوگانه جزء 
مرزی ( ۱۶۹۳5۳۲ ( مصمظ بنهمتوتمع؟۳ آهیار[ 
۵ 167160 ) مورد بررسی قرار داد. او از اینن 
روش برای الگوسازی مخزن و پی بهره گرفت و الگوی 
سد را با جزء محدود تشکیل داد [10]. ٩۹02765‏ و 
07 در سال ۲۰۰۱ اندرکنش خاک -سیال- سازه را 
با روش جزء مرزی مورد بررسی فرار دادند [11] در 
سال ۲۰۰۷ صع۷1 عصباه۷ و ۳۳۵۷۵۵۲ اثر گذاری هندسه 
منز بر بویا شاسانه اک اس ال نت زیر 
قلمرو بسامد مطالعه کردند و روش المان مرزی را در 
تحلیل به کار بردند [12]. در سال ۲۰۰۸ 0675/6006786 
و ۷1 با کمک از ووشن اخزای:سعدود کسعرش, يافعه 
20۳۳2[0) به بررسی مسأله اندرکنش سیال- سازه 
پرداختند [13]. ۲1008 و ۷2۲۵ با روش تفاضل محدود 
فشار هیدرودینامیک غیر خطی وارد بر سد صلب را 
تعیین کردند [14]. ع11 و 67 این روش رابرای 
تعیین فشار هیدرودینامیک وارد بر سد انعطاف‌پذیر 
گسترش دادند [15]. ععض با فرض ند صلب و کاربرد 
روش تفاضل محدود. مقدار فشار هیدروپویای را برای 
مخزن با شکل دلخواه محاسبه کرد [16]. 

برای تحلیل اولیه و نیز وارسی نتیجه‌ه ااز 
روش‌های ساده شده وابسته به حل بسته استفاده می‌شود. 
در سال (۱۹۹۱-۱۹۹۰) :192 و 16 یک روش نیمه 
تحلیلی برای حل مسأله اندرکنش سد و مخزن در حالت 
دوبعدی و سه بعدی در قلمرو زمان ارایه کردند. در اینن 
روش مرز دوردست در بی‌نهایت فرض شده و فشار 
هیدرودینامیک در مخزن از معادله فشار موج پیروی 
می‌کند. سد به صورت تیر اویلر- برنولی با مقطع ثابت 
قرضوی هیر انطهس زق ارانه فلع این اساتی اسفواد 
است [17,18]. از تحلیل سازه سد. ویژگی‌های پویای 


سال بیست و دوم. شماره دی ۱۳۹۰ 


ساناز محمودپور- رضا عطارنژاد 


سد. بسامدهای طبیعی و شکل مودهای ارتعاشی, برای 
کارپرد در حل مسأله اندر کنش به دست می‌آیند. از 
آن‌جایی که سازه سد با مقطع ثابت فرض شده است» حل 
معادله دیفرانسیل حرکت. به‌سادگی و از روش جداسازی 
متغیرها صورت می‌پلیرد. در عمل. سدها دارای 
مقطع‌های متغیر در ارتفاع می‌باشند. به همین دلیل برای 
موردهای کاربردی, الگوسازی سد باید به شکل تیر با 
مقطع متغیرانجام پذیرد که در این حالت تحلیل تیر و حل 
معادله دیفرانسیل حاکم بر حرکت آن نسبت په حالت 
مقطع ثابت پیچیده‌تر می‌باشد. 

فرساد و عطارنژاد در سال ۲۰۰۳ با گسترش روش 
تسای [18]» روشی نیمه تحلیلی برای بررسی مسأله 
اندرکنش در حالتی که ضخامت سد متغیر باشد ارایه 
دادند [۱۹]: در ال ۲۸۰۹ ؛ رخستی. و عظارئواد به تبحلیلن 
اندرکنش سد و مخزن در قلمرو زمان به صورت نیمه 
تحلیلی با روش تابع اولیه‌گیری حلقوی منفرد مجزا 
(056) پرداختند. سازه سد با یک تیر طره با مقطم ثاست 
و مخزن به صورت یک محیط همگن غیر چسبنده نیمه 
بی نهایت الگوسازی شد [۲۰]. 

در مقاله حاضر توزیع فشار هیدرودینامیک و حل 
معادله موج با توجه به روش ارایه شده در مرجع [18] در 
نظر گرفته شده است و روش نوینی برای تحلیل سازه 
سد ارایه خواهد شد. در این رویکرد با درنظرگرفتن سازه 
سد به صورت تیر با مقطع متغیر الگوسازی انجام می‌شود. 
روش نوین پیشنهادی برای تحلیل تیرها با مقطع متغیر 
ارایه بر پایه کاربرد تابع‌های نوینی به نام "تابع‌های شکل 
مبنا ۱۳۴9((۲ظ) مصمتاعما؟ تصمصهع‌هام‌وزا متفظ ) 
می‌باشد. اين تابع‌ها از حل معادله دیفرانسیل حرکت 
تیراویلر - برنولی با مقطع متغیر به دست می آیند. با این 
تابع‌ها ویزگی‌های پویای (بسامد ها و مودهای ارتعاشی) 
تیرها با مقطع متغیر محاسبه می‌شوند. سپس با استفاده از 
بسامدها و شکل مودهای ارتعاشی محاسباتی و انجام 
اصلاحات لازم در رابطه‌سازی مرجع [18], مسأله 
اندرکنش سد و مخزن در قلمرو زمان بررسی می شود. 

در ادامه ابتدا به شرح معادله دیفرانسیل حاکم بر 


حرکت سامانه سد و مخزن پرداخته خواهد شد و درباره 
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فشار هیدرودینامیک وارد بر بدنه سد و روابط حاکم 
توضیح کوتاهی داده می‌شود. یادآور می‌شود که این 
رابطه‌ها و معادله‌ها به شکل کامل در مرجع [18] آمده‌اند 
و در اینجا معرفی کوتاهی از آن‌ها به نظر خوانندگان 
می‌رسد. سپس روش جدید برای تحلیل سازه سد. با 
تعریف تابع‌های شکل مبنا و چگونگی به‌کارگیری این 
تابع‌هاء بیان خواهدشد. در پایان نیز تحلیل سازه سد با 
یک نمونه عددی شرح داده می‌شود و به دنبال آن تحلیل 
اندرکنش سد و مخزن با نمونه‌های دیگری انجام 


می‌پذیرد. 


معادله‌های حرکت 
با توجه به شکل (۱) معادله حرکت سد در اثر شتاب 
ناشی از زلزله و فشار هیدرودینامیک مخزن به صورت 


یر تفه یی [19 1 
2 2 2 
)۱( (۳)0,2,۸- پفلنه بس|؟ 3 ِ 


که در آن: 

مدول الاستیسیته مصالح بدنه سد 

ممان اینرسی سطح مقطع حول محور خمش (9) 
0 جرم در واحد طول سازه 

تغییر مکان سازه در راستای محور 

لا شتاب زمین در امتداد بالادست به پایین دست 
۴ فشار هیدرودینامیک 

ف بات 


1 
۳ 


3 
(۲۱۵۰ 


/ ۴3 


شکل ۱ سامانه مخزن و سد یا مقطع متغیر 


با استفاده,از روش جذاسازی سغیرهاه تفینر مکان 
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به صورت زیر فرض می‌شود: 


۳62, < _,)2(۷,۲( (۳) 
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که 9 شکل مود ارتعاش بوده و از حل معادله 
دیفرانسیل حاکم بر تير اویلر- برنولی به دست می‌آید و 
۷ مختصات کلی تابع زمان می‌باشد. لازم به یادآوری 


تسامد سامانه سازه- سیال چشم‌پوشی می‌شود. معادله (۱) 


به صورت زیر درمی آید: 
)۳( 2یا (),۳-(),۷-- (۱) ,۷ با( +() بل ,۱۷ 


با حل معادله (۳) و محاسبه ۰۷ تغییرمکان کلی از 
انطه ( انمض آهته خی 
زر شا میهف کر وا 


۱ 
( ۱۳/۹ <- ,3 
0 
نیروهای ناشی از زلزله» م۷ برابرند با: 
ط 
)0( (2) ,۵ (2) ۱0۵ [ () نا ۷ 
0 
و نیروی تعمیم‌یافته ناشی از فشار برابر است با: 


1" 
0 ,۳02 ره [ حدم 
0 


در رابطه (1) از حل معادله فشار موج به دست 
می‌آید که در ادامه بیان خواهد شد. 

فرض می‌شود فشار هیدرودینامیک در قلمرو سیال 
از معادله فشار موج پیروی کند: 
)۷( 9 - (1 ,۷۳,2 

که در آن.۳ فشار هیدرو دینامیک بدون گنجاندن 
اثر هیدرواستاتیک و ) سرعت صوت در آب می باشند. 
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دور می‌شوند. 

شرایط مرزی و اولیه در ناحیه تماسی سازه و سیال 
به این صورت هستند: 
۱. در ناحیه تماسی سازه و سیال: 


1 ,هه ۳ -[۵+)یت]م- ۳ 
91 ۳ 
)۸( 
1 در کف مخزن: 
)0( ی( 
ام 0 

۳ در سطح آزاد مخزن با صرف نظر از اثر موج‌های 
سطحی: 
۱۰ 0ب[ 


۶ با فرض اینکه مخزن در زمان صفر در حالت آزاد قرار 


داشته باشد: 

)۱۱ 0 و( 

ار 

)۱۲( 6 ۲۳۱ 
ما20 


با کمک تبدیل لاپلاس و واردکردن شرایط مرزی» 
توزیع فشار هیدرودینامیک از حل معادله (۷) محاسبه 
شده و می‌تواند به دو بخش حرکت جسم صلب سد و 
تغییر مکان سد جداسازی شود. بنابراین» رابطه (۱۳) 


نتیجه خواهد شد: 


۳, ) - ۳0+ ۳: )( ۱۳ 


0 


+ 9۵ 
(۱( 2-۲ ] ول() نا / ۶ست روت < () ۳0 
0 ا 


۳ .2 200 
کینوت ]مه 2 - )2 
0 


احط احع1 


(۱) 


) سرعت 2 در آب و ول تابع بسل می‌باشد. در 
هی هو و رات صوزخات 


دست می‌آیند: 


سال بیست و دوم. شماره دی ۱۳۹۰ 


ساناز محمودپور- رضا عطارنژاد 


/ 
02 (2) ۳ < ,04 
0 
)۱2 21-1(7) ۳ ۳ 
20 
با فرض تغییرات خطی شتاب زمین در دو گام پی 
درپی زمانی و با فرض شتاب میانگین ثابت برای سازه در 
پایان معادله (۳) به صورت زیر نوشته می‌شود: 
۲ 


2(07) ,۵ (100)2 |( یفن - (۱) ,۵2,۷ +(۱) ,۷ ,۱۸ 
0 


2 1+1 1-) 2 400 
2,٩ - 07‏ ] ات و 
0 ّ 
220 
عل[0- 6( مس مره 9 
5 1« 12۱ 
)۱۳( 
0( 5 یحظ مب حرقظ.. ررلظ ۱ 
5 0( ۷ مب ومل1 11121 
9 
(1,)۱ ,۲ ۳ 
(ا)ر] ِ )۷ ورگ 
تا نیال سک , 
)۱۸ 
که در آن: 
رو ۶ زر ۷۷ < رر۳ 
۳ ۷۷ + 1۷ < ررحط 
۳ 
3 ۵۷ < ربکا 
0-۳0۵0 ۳۳ -() :2-۷ () ر1 


همان گونه که ملاحظه می‌شود. محاسبه مقدارهای 
بالا بستگی به مقادیر ۵ و دارد که در حالت تیر با مقطع 
ثابت به صورت دقیق در دسترس می‌باشند. در حالت سد 
با مقطع متغیر محاسبه مقادیر ۵ و9 پیچیده بوده و حل 
بسته‌ای برای آن ارایه نشده است. افزون بر موردهای 
یادشده در حالت مقطع ثابت جرم نیز دارای مقدار ثابتی 
است و در روابط (۶ و۵) به‌سادگی از تابع اولیه خارج 


می‌شود. درصورتی که در تیر با مقطع متغین جرم تابعی 


سال بیست و دوم شماره دو ۱۳۹۰ 


۳ 


از ارتفاع می باشد که به پیدایش تابع اولیه‌های پیچیده‌تر 
می‌انجامد. در بخش بعد محاسبه مقدارهای بسامد طبیعی 
و شکل مودهای ارتعاشی برای تیر با مقطع متغیر شرح 
داده خواهد شد وسپس سد به صورت تیر با مقطع متغیر 
الگوسازی و تحلیل می‌شود. در مقاله حاضر در کنار 
روش نوینی که برای تحلیل پویای تیر با مقطع متغیر ارایه 
شده اصلاحاتی نیز در حل معادله (۱۷) انجام شده است. 
این اصطلاحات بر پایه کاربرد دنباله‌های توانی و بسط 
تابع‌های چندجمله‌ای برای تابع اولیه‌گیری تابع‌های بسل 
در رابطه (۱۷) هستند. بنا به رابطه (۱۷) تابع اولیه‌گیری از 
تابع‌های بسل ظاهر می‌شود. در مرجع [18] از تابع‌های 
۵ به جای تابع اولیه‌های تابع بسل استفاده شده 
است و پاسخ‌ها با ۳۵ جمله در معادله (۱۷) محاسبه 
شده‌اند. در روش اصلاحی پیشنهادی با ۱۰ جمله پاسخ‌ها 
با دقت خوبی محاسبه می‌شود. این اصلاح در روش حل 
دقت و سرعت را بالا می‌برد. 


محامبه شکل مودها و بسامدهای تیر با مقطع متغیر 
با توجه به اهمیت محاسبه شکل مودها و بسامد تیر با 
مقطع متغیر در اين مقاله روش نوینی برای محاسبه آن‌ها 
ارایه می‌شود. به اين منظور لازم است تا ابتدا به معرفی 
یهاش کل مسا پر واه شود اتفریت‌های م ری یه 
تابم‌های شکل مبنا که پایه و اساس روش پیشنهادی است 
۵ ون )سل اند 

با توجه به اصل تقابل کار بتی- ماکسول. می‌توان 
تابع‌های شکل مبنا را مانند جدول (۲) تعریف کرد. 

با توجه به این تعریف‌ها ماتریس نرمی گرهی با 
رابطه‌های (۲۳) و (۲۶) به دست می‌آید: 


| (00,ظ‎  0,)0( 

۲۳۱( بط ارطه اع رظ 
ما که مب 41 ۱ 
بیط ل(لایط | 

(۲( روط م0 | مر 
رآ 1 3 1 ۳ 
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ماتریس سختی گره‌ها نیز به این صورت بیان می‌شود: 


(۲۵( ۳ رد ۱ 0 5-ه 
میک 0 


با معرفی 9 به صورت: 


۲0( و عبط با بطاعه 


۷1 

و بهره گرفتن از مفهوم‌های روش اجزای محدود. 

ماتریس‌های جرم و سختی سازه‌ای با رابطه‌های (۲۷ و 
(۲۸) محاسبه خواهند شد: [21] 


ن1 
۳( ی [ ٩‏ < ۱ 
0 


۱۳ 
۳۸( هر [ 1 1 
0 


با توجه به ماتریس‌های جرم و سختی سازه تحلیل 
سازه و تعیین دیگر عامل‌ها به‌ساد گی امکان‌پذیر خواهد 
در حالت ایستا برای تعبین جزء‌ی بردار 9 که در واقع 
نتیجه حل معادله دیفرانسیل تیر در حالت ایستا است. 


بهره گرفتت, 
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۱ 1 تعریف نشانه 
| و تخیر مکانقائم گره سمت چپ تیر طره در اثر 
5 گ ۱ بار واحدی که در فاصله اثر می کند . 9 
1 ۷ 
گر | زاویه چرخحش گره سمت چپ تير طره در اثر 
1 2 بار واحدی که در فاصله اثر می کند. 9 
9 تغییر مکان قائم گره سمت راست تیر طره در اثر 
2 
۱ 1 2 بار واحدی که در فاصله اثر می کند. 9 
۷ 
ظ زاویه چرخش گره سمت چپ تیر طره در اثر 
02 ِ 
1 کر بار واحدی که در فاصله اثر می کند. و 


اند رکنشر سا و محزن در قلمرو زمان 


در حالت پویا باید معادله دیفرانسیل حاکم بر 
حرکت. رابطه(۲۹) را حل کرد. 


22۷7۷6 ,62۷7۷6 7 
۳ 06(۰) ۵ + ۱ 7 مت 20 


(ا)0 < ِ 


(۲۹( 
با فوض: 
(۳۲۰) (اه)ومنه + را هصرع (00 ع ()یع06 


(۳۱ (10 )مه + (6818)0۵1)() ۷ < (ای) ۷ 


معادله دیفرانسیل حرکت به شکل زير نوشته می‌شود: 


۳ 
ما 


170 
۳ 


) ۵.0۸00۰ 0 ۳۲( 


می‌توان روش‌های عددی گوناگون را برای حل این 
معادله دیفرانسیل به کار برد. پس از محاسبه جزء‌ی بردار 
9و تشکیل ماتریس‌های جرم و سختی سازه می‌توان 
ریاد مدای نی با منم موز تعاس 
کرد. در اين مقاله از روش تبدیل دیفرانسیل برای حل 
معادله حرکت و محاسبه اجزای بردار 9 استفاده می‌شود. 
با ون تال زب شک رز تقایل دشراشنیز راز 


(۳۲) به صورت زیر نوشته خواهد شد: 


سال بیست و دوم شماره دو, ۱۳۹۰ 
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۳0 + 1()+ 2(6+ 3()6 + 4( ۲۷۷ + 4( + 


1-0 


290-3 + 1(۳۲10 -1+1(0+ 10+ 2(6+3(۲۷۵+3(+ 


0ع[ 


(1(6+2(۷۷۵+2 + 2(6 + 2(۳۴10 +1- 0( جوز ۱ 


- 2۸0 -1(۷۷( 


(۳۳ 


با از رابطه (۳۳) همه جمله‌های ۱۷6۵ به ازای 
4 بو خسن ههار صماتتنه ازل::(1۷)0: 
(۰1۷0 (۷۷62 و (۷۷)3 محاسبه می‌ شوند. بنابراین برای به 
دست آوردن تغییر مکان محاسبه چهار جمله اول کافی 
هست. با وارد کردن شرایط مرزی مناسب می‌توان 
جمله‌های نامبرده را محاسبه کرد. با توجه به تعریف ارایه 
شده در جدول (۲) برای تابع‌ها شکل مبنا . شرایط مرزی 


مناسب برای هر تابع شکل را می توان به صورت زیر در 


نظر گرفت. 
2 بعا۳ 0 ۳ 21 م۷ : 0 
(۳) 
2 9 وت با ند م۷ 29 # : ,و0 
(۳۵( 
سس ۷ 2 بخ 0 م0 با مب ز یروط 
۳0( 
1 ۷ 0 ۷ 20 0 0< میا :و9 
۳۷( 
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سال بیست و دوم شماره دو ۱۳۹۰ 


جدول ۲ تعریف تابع‌ها شکل مب با توجه به اصل تقابل کار بتی-ماکسول 


تغییر مکان قائم در فاصلهر تحت تاثیر بار واحدی که در 
گره سمت چپ تیر طره وارد می شود . ۷1 

7 زاویه چرخش در فاصله تحت تاثیر لنگر واحدی که در 

: ۱ ۱ 1 

و9 ۹ گره سمت چپ تیر طره وارد می شود . 
تغییر مکان قائم در فاصله تحت تاثیر بار واحدی که در 
گره سمت راست تیر وارد می شود . ۷2 
زاویه چرخش در فاصله تحت تاثیر لنگر واحدی که در 


گره سمت راست تیر وارد می شود . 


4۵ 


با کاربرد تئوری تير اویلر- برنولی زاویه چرخحش» 


محاسبه کرد: 
"1 
۳۸0( 73 000 
(۳۵) 3 < ()/۱ 
1 
9 4 
(۶۰) ۱ ت ی < (۶) ۷ 


با این رابطه‌ها؛ چهار جمله اول 1۷66 و در نتیجه 
تغییر مکان تیر محاسبه می‌شود. 

در ادامه برای دریافت بهتری از چگونگی کاربرد 
تابع‌های شکل مبنا در محاسبه ویژگی‌های سامانه یک 
نمونه عددی برای محاسبه بسامد طبیعی و شکل مودهای 
ارتعاشی تیر با مقطع متغیر ارایه خواهد شد. در نمونه ۲ 
برای وارسی درستی صحت روش پیشنهادی. محاسبات 
مربوط به سد با مقطع ثابت با روش جدید انجام شده و 
نتیجه‌های حاصل با نتیجه‌های موجود در نوشته‌های فنی 
مقایسه شده اند. در نمونه ۳ نتیجه‌های تحلیل پویای سد 
و مخزن در حالت تير با مقطع متغیر انجام شده است. 
نتیجه‌های حاصل با نتیجه‌های حاصل از برنامه 8۸0۲ 
که بر مبنای روش اجزای محدود است [6 روش نیمه 
تحلیلی [۱۹] و روش تفاضل محدود [۲۲] مورد بررسی 
و مقایسه قرارگرفته‌اند. بررسی نتیجه‌ها بیانگر کارایی و 
درستی روش پیشنهادی می‌باشد. 
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نمونه ها 

نمونه 1 
در این نمونه چگونگی محاسبه تابع‌های شکل مبنا برای 
یک تير با مقطع متغیر با روش تبدیل دیفرانسیل شرح 
داده شده و ماتریس‌های سختی و جرم. محاسبه خواهند 
شد. نتیجه‌های حاصل با نتیجه‌های موجود در نوشته‌های 
فنی مقایسه و فرضیات زير در نظر گرفته شده‌اند: 

4 > بآ ز ظ / 2700182 م زک / 71*101 - ظ 

۸) 2۱2۱0 0-058( 

سطح مقطع و لنگر لختی مقطع تیر به صورت زیر 

تغییر می‌کند: 


(1-0.58) 21*107 (168 و (01-0.58 1.27107 )۸ 


چهار جمله اول به این صورت تعریف می‌شود: 


(۱) ۲ - (0) 1۷ 
(6۲) 5 )۷۷ 
(۳: - (۷۷۵2 
(1ع) 7 - (3) ۷۷ 


بافرض 0 - »و جای‌گذاری رابطه‌های (41-4۱) 
در معادله (۳۳) این رابطه به دست فش اب 


(۵ع) 52 + 0.125 - 0.1947 - (4) ۱۷ 


به همین ترتیب با فرض 21 6 و جای‌گذاری 
رابطه‌های (۶1-6۱) در معادله (۳۳) داریم. رابطه (۶1) 


حاصل خواهد شد: 

(7ع) 572 + 0.10 - 0.03895 + 0.15587 - (5) 1۷ 

به ازای 2 < 16: 

۱۷۷ )6( - 0.09747 + 0.03259 -0.04950: +0252  )۷( 
به همین روش می‌توان جمله‌های بعدی را نیز‎ 


محاسبه کرد. در این نمونه ۱۸ جمله در نظر گرفته 
می‌شود. با کمک ویژگی‌های روش تبدیل دیفرانسیل تابع 
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تغییر مکان به صورت زیر نوشته می‌شود: 
(6۸) () 2,۷۷ - )۷ 


با وارد کردن شرایط مرزی برای تعیین ,9 این 
چهار معادله را می‌توان به دست آورد: 


2130010 + 426007 <2 4 
۷-0 
1.5707۲ +1.11925 + 0.67460 + 2.75737 20 
3.1463۲ +1.78275 + 0.29230 + 12.260372 20 


39 
با حل دستگاه معادله .4٩(‏ مقدارهای 2 5 دا و 7 
25 0-20 0.0220 <5 0.0183 < ۲ 
(6۰( 


در پایان ,, به صورت زير محاسبه خواهد شد: 


9. --1.834562+2.2044-726 + 1.50237 
2 4450-3٩-1 3229-367 -6.9495 4۲ 
-3 5108647 -1.83696-421-9.70586-5 
51259561-2 69535 *عکع1.4101-‎ 
-7 34076-627 2.80486 -1 9645-67 
- 1.01096 51871-767 265487۲ 

5۱( 


به همین ترتیب» دیگر تابع‌های شکل مبنا نیز 


محاسبه خواهند شد: 


7ْ150-3-5.66796-26-1.12686-3 5.51 < وط 

*ع 496934 + *4-ع9142. 7+ 5.6338047- 
*ع7.23286-5+ 7ع1.37480-4+*ع2.6758-4+ 

اع 1.0675 +2.02766-5+ *ع5-ع8326. 3+ 
*ع2.91176-667+1.50996-6+ 58930-617 5+ 
0607-۳7 2+ 01496-761 4+ 3 79947 7+ 
*ع7-ع1.0549+ 

(۲( 


سال بیست و دوم. شماره دی ۱۳۹۰ 
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47ع1 8444+ 367ع7.0670 2 بط 
"388504 3.26716-47-1- 2.50064 
1.06340-5- 2.36415 5.59770-567- 
15286-6۳7 2.37746-61-1- 4971506-66۳1 
3737-77 2-1 2.77786-7- 64150-767 5- 
*6.81420-866-3.38786-867-1.68716-6- 
(5۳( 
6.07646-4+ *2.08150-36 < روط 
90 7ع6.52920-5+ *ع9555-4. 1+ 
*-20340. 1+ 16-58 2.077+ 35810-567. ۵+ ؟ع5- 
97106-62۳ 3.66906-61+1+ 1ع6.72176-6+ 
8781-77 2+ 50767 5+ 04666-627. 1+ 
*-3.98576+ 7.73536-8617-+ 49550-72۳ 1+ 
(( 


ماتریس‌های نرمی گرهی با رابطه‌های (۱۸) و )۱٩(‏ 


ند تشکیا می‌شوند: 


- 5 
11 1 ۱ )6۵۵( 


0.0055 4 


"(2 (۲ 1 
22 - ۱ (51( 


0.0031 5 


2024 79 

11 < (۷( 
1187. 5 

و 79- 5678 ] 
دور 9 روت | 2۳ 


با استفاده از رابطه (۲۵) می توان لوششت: 


302.4 1187.9 ۰.0 0 

11879 39525 0 0 

1 ۵ 0 567.8 -79 
0 0 9 7 


(5٩( 
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۷ 


با توجه به رابطه (۵۹) و به کار بردن رابطه‌های 
(۲۸-۲۷) ماتریس‌های سختی و جرم محاسبه می‌شوند و 
با کمک رابطه (1۰) بسامدها به دست می‌آیند. در پایان 
شکل مودها محاسبه خواهند شد. 
(.) ۷>۵. ۳-۵۶ 

در جدول (۳) بسامد طبیعی با روش مورد بحث 
محاسبه و نتیجه‌های حاصل با نتیجه‌های موجود در 
مرجع [23] مقایسه شده‌اند. 


جدول ۳ بسامد بدون بعد تیر با مقطعی که به صورت خطی تغییر 


می کند 
مرجع [23] | روش حاضر | بسامد بی بعد 
9 3.82378 مود اول 
3 18.31726 مود دوم 
8 | 47.26482 | مود سوم 
05 ۳ 9045047 | مود چهارم 
2 ۳ 148.00174 | مود پنجم 


نمونه ۲ 


در این نمونه پاسخ پویای سد با مقطع ثابت با تابع‌های 
شکل پیشنهادی محاسبه و نتیجه‌ها با حل تحلیلی 1921 
6 مقایسه شده‌اند[18]. ضخامت سد ثابت و برابر 
۲ ارتفاع سد ۱۸۰۲0 جرم واحد طول آن ع/ع۳/۱۵0 
و *۲ ,۹/۸:۳۷۱۵۲ عآت فرض شده است. 

در جدول (۶) بسامدهای طبیعی تیر با مقطع ثابت 
با کمک روش پیشنهادی محاسبه شده و نتیجه‌های 
حاصل با نتیجه‌های دقیق مرجع‌های [23] و[24] مورد 
بررسی و مقایسه قرارگرفته‌اند. در شکل (۲) پاسخ پویای 
سامانه سد و مخزن در برابر زلزله السنترو که در مرجع 
[18] ارایه شده, آمده است و نتیجه‌های حاصل از تحلیل 
با استفاده از تابع‌های شکل و مقایسه آن‌ها با حل تحلیلی 
1 در شکل‌های )٩-۵(‏ مشاهده می‌شود . 


نشریه مهندسی عمرال دانشگاه فردوسی مشهد 


7۸ 


۱۱۱۹۹ 


و عم وا و۲ 


اند رکنش سا و محزن در قلمرو زمان 


۱۱ 


تنشریه مهندسی عمران دانشگاه فردوسی مشهد 


شکل ۳ سیستم مخزن و سد با مقطع ابت 


111۳066660 


شکل ۶ شتاب رمپ 


0.02 6 


‌ 


۹ 


سال بیست و دوم. شماره دی ۱۳۹۰ 


ساناز محمودپور- رضا عطارنژاد 


مرجع [24] 
2-1602 
2.5( 
6۱972 
120902 


19986 


جدول ۶ بسامد بدون بعد تیر با مقطع ثابت 


مرجع [23] 
2-1602 
7-25 
6۰972 
120902 


19986 


]1۳08)586( 


روش حاضر 
2-101 
2-2-0( 
2۱-2-91 
120۰-90191 
12-3-9334 


شکل ۵ مقایسه تغییر مکان تاج سد با مخزن پر تحت اثر زلزله ال سنترو 


شکل 7 مقایسه تغییر مکان تاج سد با مخزن خالی تحت اثر زلزله ال سنترو 


سألم یب بیست و دوم شماره دو ۱۳۹۰ 


]1۳08)580( 


801/۵ عععاموزط 


8/۵ صعوامون0 


14 


نشریه مهندسی عمرال دانشگاه فردوسی مشهد 


تشر 


ب 


مهندسی عمران دانشگاه فردوسی مشهد 


سبال 


بیست 


و دوم. شماره دی ۱۳۹۰ 


۲1۳۱۵)686( 


۱ 
3 


۱ 
۳۹ 
9 
ِ‌ 


۱ 
‌ِ 
اح 
ب‌ِ 


9/0 ساوناج(سمع) 


ی :۳۱۳ 


شکل ۸ مقاٍ 


۸ مقایسه 


ف یر و 


71 


ثر 


شتاب رمپ 


]1۳8)8860( 


03 


ی زز ( ۲ 


28۱0و وناز زنلاع) 


00۴5 سس 


و تم 0.05 


شکل ۷ مایسه تغییر 


تحت ار 


بر 


زلرله ال 


نشرو 


(0وو)ع1۳8] 


3 


دنس 9و 


۱ 
۹1 
۵ 
ب‌ِ 


0.01 


نان 


۳ 


0و 9دوساونج(نتم6) 


002 


۱۵۵9۳۷۵۱-۵۵۴ ام 
5 ۲95۶۳۷۵۱-8۵۴ الا 


003 


اندرکنش سد و مخزن در قلمرو زمان 


ساناز محمودپور- رضا عطارنژاد 


نمونه ۳ 

در این نمونه پاسخ پویای سد با مقطعی که به صورت 
خطی تغییر می‌کند (شکل ۱» با تابع‌های شکل ارایه شده 
محاسبه و نتیجه‌ها با نتیجه‌های حاصل از برنامه 3۸0 
[6] و روش تفاضل محدود [۱۹] و روش نیمه تحلیلی 
[۲۲] مقایسه شده‌اند (شکل‌های ۱۰- ۱۵). ضخامت سد 
در کف مخزن "۲۵۲ و در بالای سد ۱۲/۵۲9۲ ارتفاع 
سد ۱۸۰۲ چگالی آن ۱00/۲ع/۲ و ضریب کشسانی 
۳/۴ فرض شده است. 

نمونه ۶ 

در اين نمونه اثر در نظر گرفتن تغییرات مقطع در پاسخ 


]1۳08)6960( 


۷ 


سامانه سد و مخزن بررسی شده است. دو حالت مخزن 
با مقطم ثابست و ضخامت ۲۵۲۵ در کف و مخزن با 
ضخامت متغیر (ضخامت ۲۵۲ در کف و ۱۲/۵۲ در 
بالای سد) درنظر گرفته شده و پاسخ پویای سد و تغییبرات 
فشار هیدرودینامیک در دو حالت با یکدیگر مورد مقایسه و 
بررسی قرار گرفته‌اند (شکل‌های (۱۳) و (۱۷)). 
همان‌طورکه ملاحظه می‌شود. درنظرگرفتن تغییرات 
مقطع. تاثیر تعبین‌کننده‌ای در پاسخ سامانه دارد و الگوکردن 
سد با مقطع ثاببت درجهت اطمینان نیست و سفارش 


005 
اوعااراوطه اعد 
بت 
سم | ومع ددعت [008 


0۲۵ ححه 


(و/) اطع ععواموزط 
[- 


شکل ۱۰ مقایسه تغییر مکان تاج سد با مقطع متغیر در حالت مخزن پر تحت اثر زلزله ال سنترو 


]1۳8)9860( 


اوعا 52۳0-202 


ارت ۱0 


شکل ۱۱ مقایسه تغییر مکان تاج سد با مقطع متغیر در حالت مخزن خالی تحت اثر زلزله ال سنترو 


سال بیست و دوم شماره دو ۱۳۹۰ 


نشریه مهندسی عمرال دانشگاه فردوسی مشهد 
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شکل ۱۲ مقایسه 
شکل ۱۳ مقایسه 


مهندسی عمران دانشگاه فردوسی مشهد 


اسست ۲۶ 
اسف ۱۳ 


شکل ۱۶ مقایسه تغییر مکان تاج سد با مقطع 


]1۳08)586( 
]1۳08)886( 
]1۳8)6860( 
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۱ 
بط 
بِ 


۱0 ونر زننمع) 


اجع1 52۳0-020 


۷۲ 


اند رکنش سد و محزن در قلمرو زمال 


سبال 


بیست 


و دوم شماره دو ۱۳۹۰ 


نشر 


2 


مهندسی عمران دانشگاه فردوسی مشهد 


شکل ۱۷ مقایسه تغییرات فشار هیدری 


دینا 


میک سد با مقطع 


تا 


بت و 


از 


تحت از 


زلزله السترو 


]1۳08)90( 
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ود عوه۲ع واطو۷۲ 


و2۲۵ اطاحاوده6 


شکل ۱۱ مقایسه تغییر مکان تاج سد با مقطع ثاب 


بت و 


اثر زلزله السنترو 
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ات جات ع ای سدع 2 17 
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0 2 فساوناچ(ن6) 


طمااعوو عوه۲ه وااج[۱۷2۱ 
۴ 6۲۵95 060۳091301 


۵ مقایسه 


مکان تاج سد با مقطع متغیر 


در حالت مخزن پر و خالی 


ه رد 


7 


بر 


شتاب رمپ 


]1۳8)690( 


و 


ماه ‏ عت هت ات ۱۴ 
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هت با تعاس میم مره له مه سا 
۱۳ 


1 
0۳ 
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واولهرقدوسیاویجزنن8۸) 


ساناز محمودیور- رضا عطارنژاد 


۷۳ 


۷ 


نتیجه گیری 
در این مقاله تحلیل دقیق پویای اندرکنش سد و مخزن در 
قلمرو زمان بررسی شد. سازه سد انعطاف‌پذیر و مخزن تا 
بی‌نهایت در نظر گرفته شد. الگوسازی سازه سد به 
صورت تیر اویلر- برنولی با مقطع متغیر و درنظرگرفتن 
بدنه سد در بالادست به صورت قائم انجام پذیرفت. 
روش نوینی برای تحلیل تیر با مقطع متغیر ارایه شد. این 
روش نوین بر مبنای کاربرد تابع‌های جدیدی است که 
تابع‌های شکل مبنا نام دارند.این تابع‌ها از حل معادله 
دیفرانسیل حرکت تیراویلر - برنولی با مقطع متغیر به 
دست می‌آیند. با استفاده از اين تابع‌ها ویژگی‌های پوبا 
(بسامدها و مودهای ارتعاشی) تیرها با مقطع متغیر 
محاسبه می‌شوند. سپس با کمک بسامدها وشکل مودهای 
ارتعاشی محاسباتی مسأله اندرکنش سد و مخزن در 
قلمرو زمان مورد بررسی قرار گرفت. در فرآیند حل 
مسأله اندرکنش» اصلاحاتی صورت گرفت. در این 
اصلاحات دنباله توانی و بسط تابع‌های چند جمله‌ای 


برای تابع اولیه گیری تابع‌های بسل به کار رفت و با 


و وت 


اند رکنشر سا و محزن در قلمرو زمان 


تعداد جمله‌های کمتر جواب‌های دقیقی شاه دسعت. امد 
در حالت مخزن خالی بررسی ماو ها ان بسا 
نتیجه‌های موجود در نوشته‌های فنی مقایسه شد. مقایسه 
نتیجه‌ها پیانگر کارا و درستی روش پیشنهادی 

در ادامی اندرکنش سد و مخزن در قلمرو زمان و 
در اثر زلزله ال سنترو و شتاب رمپ مورد بررسی قرار 
گرفت. شایان توجه است که نتیجه‌های تحلیل اندرکنش 
سازه و مخزن, افزون بر آثر حضور مخزن به تغییر 
شکل‌های سامانه سازه وایسته می‌باشد و از اثر نغییر در 
مقطع سد نمی‌توان چشم پوشید. از این رو روش دقیق 
تعیین تغییر شکل‌های سامانه سازه در حالت مقطع متغیسر 
ارایه شده است که این امر امکان تحلیل اندرکنش سازه- 
سیال را با دقت بالا فراهم می‌سازد. در پایان نتیجه‌ها با 
نتیجه‌های موجود در نوشته‌های فنی مقایسه شده‌اند که 
این مقایسه. دقفت بالای محاسبات و کارایی روش 


پیشنهادی را نشان می‌دهد. 


1933(۰) ,418-433 .00 ,98 رفظ .عم ماه مصسط عم من فمتاعوم0۲ ۷۷۵۷۵ ,۲۰۸۲ رلت2هعماوه ۷۷ . .1 


وتگن) کر ۷۵6( .۳909 ]۵ .[ .کم هن اجه عصتسل وق وم فعتاعوعم متصهوصررل۳۱/0۲۵0 تونک به رقت0مط) ...2 


93)6(, 00. 205-223, )1967(, 


0۳ ۳:۵۵ .2 ۳۵۲۲ .وماهامطانهه عصتسل هل مصامماه وم وماووهيم متصصهه و0 رنه هن .3 


,(1978) ,343-348 .وج ,87 ,۱۸۵6۲ با ۵۲ .ل 


صنط وه ۵۶ جمناه:۷10 معط ۶۵۲ فصمتبتامو کصعجصماه -0ها بو جرج )مقر تر.ظ رعصم]- همه ,خه۷ بق۲۵ وی بلع]۷ .4 


.(79-88,)1979 .00 و(1)2 و.عع؟ طوعع0 ,۸۳۴۳۱ 11001۰۳ 6و فطه م 56260 ممبتامتاو متاقوعاه 


۰ ,165 .۷۲66 0زیاا۳ .ل ".عملهنمطصانوه عصس فصهل م11 مم فمتناوعمزن متمصوصرول۳۱۷0۵ .نا مینالنا. .5 


131-145, )1986( 


۳ 2۳21 مامتمط۵ه 0 ولوامصج ملقنامطانهه ۶۵۲ هنومن تمانام‌جمن ۳۸ ویک نهر 00۵18 ضرع وعتصو۳ ...6 


,(1984) :۱۸۵۲-84 رومامیاهظ رهتمگتلمت ۵۶ اوهتنصته ۱۲۳۲/۳۳۲۹۵۲۵4/11 ۲۱۵۵0۵۲۲ 


تشریه مهندسی عمرال دانشگاه فردوسی مشهد 


سال بیست و دوم. شماره دی ۱۳۹۰ 


ساناز محمودپور- رضا عطارنژاد 


موی مامتان ۲ 12-1۵501۷01 گم وتورافمم تاماهتا ملقتامطنهفط ر.ظ.۲ پبارلظ مه ی فطل 
,3 .5۱99 عامبواع6 ۵ برجم آمصمم۲ .۲.6 کف ]۵ ون0ععع0ظر ۱۵۱۵۵۵۰۳ اصعصعاظ مهم عصنوت 
1990(۰) ,85-94 .00 

ولو افص ممتامجماص نت۲ ۲۵۲ عصتآمنام ۳۳۱۷۲-۰۲۳۸۲ و۵" را وعاصم 20 ۷۵ 0۵۰ 0۲۶اوظ 
1991(۰) ,00۰1151-1168 و31 ,.۲۳۱۳۵0 ۷۵۷۰ ۳۳۵۱۱۵۲۰ ]۲.0 ۳۲ ۲بصتعحطمل محصت عط 

-5011 لقصمزعومصصز(2-1 ۲۵۲ ۱۷۲۵061 اصمحصما ماتصزص]-ع) ۳۲۱1 املاهت۳2 ۸ ویک مصتععک 0ص2 ب۲1۳ رلله ۷ 
2003(۰) ,249-253 .00 و23 ,او موه زر آزمی ۲ .عصصمامامتج ممتامعتماصا سامتاه 

۸ مواآمنامع ها ممتامجماص ۲ز1۵90۳۷0- و1 ۵۶ ولو رواقمم متصصهطرنا کصمتفمه1 1 رو رصهاوتم‌لنمیی۴ 
2004(۰) ,692-707 00۰ م21 و0۳۵ 8989 ۲۳۱۲۲۱ 208 

7 عصصمامامم ممتامهعاصا معتتامناتای[01و نب ۵۶ ولو امه متصهطونا و .۲۷۷۰۲ وتتاعصه]۱ 4ص2 ...1 رگ وعتوم8 
.(12,)2006 498-5 .وو, 219 روعز۳5 0۵0۳۵۷۸ ]0 .[ ,"00اعحظ کحمصرعآه بتع‌صتامها فص 

تمد معط وم امومع ۲۵۹6۲۷01۲ ۵۴ قامع]گه م1 ,۲ .۲ ]۳۲۵۷۵5 28 با ۷۰ قصتاه ۸۵,۷ ۱۷۲ مصهان۱ 
2007(۰) ,00۰1441-1459 و36 و 1۰ 50۳6 وحوح ع و۲۵۳ و" عصهل کته 0۴ ۲۵500856 

تاه معصمتم۵۵0/]۵ظ۱۷]۵۱ )صمجوماظ مانصز۳ ۱۵060 مه ۲ ب۸ عمهعگ۵۱ ۲۷۷ ,۲۷۷۵11 24 .مرتمعهطاجماوتعت0) 
-1699 .وو, 197 .۲۳9۳9 ۲۷۵۵۲۰ .امرصزظ ,۱۷/۵۵۱ یامن رت صمتاممهاصا متیام‌تاای۱10؟ ۶۵۲ طمومتمرجر2 2960 
2008(۰) ,1714 

۰ ۱۵۵00 ون ۵10 و۵ ماووهتم متهصرول0ت را تممامما را رقمه ۲۷۷ 20 کل 1 رقصتا۲ 
.(1987) ,482-499 .00 و113 .۷۲66( .۳۲۱۱90 

۰ رل 126 و ۵‏ وتا8و ۳۱۵‏ ملصوصرول م۳۱۱0 تقمصتمما. ر.۳.ظ رقف 24 کهآ ما۲ 
(1372-1391,)1990 00۰ ,116, (نکه).:/۷۲۵۵( 

۲ اوتااصانه 1 وعمه ون وم معناقوهع۳ . متصصمصرول متا تمممتامم . .۲.ظ یصفط) 
404-413,)1994(۰ .00 و32 0۲۵۵۱۱6۵560۲0 .ل رد وهم2طه 

۶ کلوفصح صتقصم-عمحصت به۶ 0مطامحظ هم /وامصه-تصطمو خر "ولیک رتاک 2040 ).0 رععنآ 200 .ق وله15 
(19902) ,913-933 .00 ب(29)5 .۳98 ۷/۵۵۰ ۱۳۵۱۱۵۲۰ ۵۲ ۱ و عطمتامحعاصا 0 0و0صصرع 

بل :5۵10100 )۲۶۵۵ :۲ ابو موه تز0 0و۵ صصهل ۵۶ وعوررلمصه صتقصرم عمط ۳ ریق رلقو] 200 ,0.0 ر8عن[ 


1990-2006,)1991(۰ .00 و(117)9 رتلبن)کظ .۷۲66 .۳990 ]0 


10 


11۰ 


12 


13. 


14 


15. 


16۰ 


17. 


18. 


٩‏ فرساد. علیرضا. و عطارنژاه رضا " حل بسته اندرکنش سد و مخزن در قلمرو زمان با در نظر گرفتن تغییرات ضخامت سد" نشریه 


دانشکده فنی» جلد ۳٩‏ شماره ۳ ۳۰-۳۲۹ (۱۳۸۶). 


۰ کریمی. کیانوش و عطارنژاده رضا " اندرکنش سد و مخزن در دامنه زمان به وسیله روش تابع اولیه حلقوی منفرد مجزا (۲560)" 


نشریه دانشکده فنی. دوره ۶۲ شماره ۶: ۵۰۱۸-۶۹۷ (۱۳۸۷). 


۰ 02۳006۰ متامصصوزرمما ۵۶ ولورافمم مه فصمتامصیا؟ اممصصمم‌هام‌وزنا متفوظ ریک ر2[صها۸ .21 


(2010) ,733-745 00۰ ,27 ولا009) 


۷۵ 


سال بیست و دوم شماره دی ۱۳۹۰ نشریه مهندسی عمرال دانشگاه فردوسی مشهد 


۲ عبذاللهی: شجمود. ۲ تیل بویا انتزکیشنسندت مشون مد مین رون یمه تخلیای بایان تاه کارشتناشی: اوشنه دانشگاه کهران؛ 
دانشکده مهندسی عمران ۱۳۸۹(۰) 

فص عصنعه عصفهه 2۵۵۲60 مصتاهاوز ۵۶ ممتاه۷۱0 م۲۲ روگنا مطمیامع مضه بط وتا ویگا.ل ر6ع206[6ظ ,23 
1034-1054,)2006(۰ .00 ,298 ر۷:۳۵80۳ ۵۳00 509۳۵0 ]0 ٩.‏ ,800:۲ظ)عظ وعمطگگلاو متصعطول 

۷ ,۱4مطصنم 4صمتوملا جد۷,صمتانکظ 3۳ رد کممتامانروی مه معط دینومرظ آم یت۱۷۲۰۲5۲۳ و22ظ :24 
(1991) ,۷0116 

-01 ۶0۲ ومم1تاقصه قققعط )صماعتعومع صرح ففمصتاد متتامصطممع مطا ۵۶ ممتاه001۷ مط) ۵ ریک رله[مصتها۸ .25 
لمصماهانامصمن ۵ فقمتمصمن) صوعم مت ۵۴ ففصتلمعع0۳0 1 و قتامم‌صصرماه حصقهط تالناممرهظ-عماباظ متقصصفاتم 
2000(۰) 2۲60108 و.عطعطط فصح ,نع ,أمورض ط1 بطاه]۱۷ 

۱3۹ 0 15۳۲ ۶ کعط۳۱۵۵۵60 ط1 رو" .عصصععظ )یرصم که صم۷۱0:2 ۲۳۲۲۵۵رب؟۳ ر20زمصهتاظ .26 
2002(۰) و۷01 6۷ مرجم ,۱۷۲۵0 


نشریه مهندسی عمران دانشگاه فردوسی مشهد سال بیست و دوم شماره دو» ۱۳۹۰ 


